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1. Introduction 


Iron and nickel nuclei are normally assumed to be the most stable form of hadronic matter at 
zero external pressure. In principle, this should be testable from the basic theory of strong inter¬ 
actions, Quantum Chromo Dynamics (QCD), but in practice this is impossible in the foreseeable 
future. QCD is not suited for finite density calculations and calculations involving many degrees of 
freedom, so much of our theoretical knowledge about dense matter is based on phenomenological 
model calculations that try to incorporate and parametrize some of the main features of the strong 
interactions, such as confinement and asymptotic freedom. In some of these studies, notably stud¬ 
ies based on the MIT bag model, it has been shown, that there is a significant range of parameters 
for which a three-flavor quark phase with roughly equal numbers of up, down and strange quarks 
(called strange quark matter in bulk, and strangelets in small lumps) has lower energy than a nu¬ 
cleus with the same baryon number. Thus, strangelets rather than nuclei may be the ground state 
of hadronic matter llT],|i||,|]. 

A range of questions immediately occur if strangelets are more stable than nuclei: 

1. Why are we here—why don’t our nuclei decay? 

2. Can strangelets be created in the laboratory? 

3. Can strange quark matter be found in space? 


The answers to all three questions rely on the properties of strangelets summarized in the 
following. The answers to the questions (some of which will be further explored below) are briefly: 


1. Very low baryon number strangelets are likely to be unstable, even if strange matter in bulk 
is absolutely stable, and for intermediate baryon numbers the transformation of a nucleus 
requires an improbable high-order weak interaction to simultaneously transform of order A 
up and down quarks into strange quarks, where A is the baryon number. This explains the 
stability of our nuclei. Strangelets have a positive electric charge, and therefore repel nuclei 
from the surface, thus minimizing the risk of growth via absorption of nuclei ^ ^ |^] . 


2 . 


In principle strangelets might be formed in ultrarelativistic heavy-ion collisions by coales¬ 
cence or through a distillation mechanism leading to strangeness enhancement. However, the 
available baryon number is low, which makes it difficult to cross the low-A stability cutoff, 
and furthermore the environment is hot, so the process has rightly been compared to the cre¬ 
ation of ice cubes in a furnace. Experiments have been performed, but with negative results 


0111’HI- 


3. The most promising place to search for strange quark matter is in the cosmos. Originally 
strange quark matter “nuggets” surviving from the cosmological quark-hadron phase transi¬ 
tion were suggested as a possible and even natural candidate for the cosmological dark matter 
1^], but later studies showed, that this was probably unlikely due to the high temperature en¬ 
vironment that would lead to evaporation of the nuggets [|^, |T^ [^, [^, |^. Recent 

lattice results indicating, that the quark-hadron transition at low chemical potential (like the 
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early Universe) is not a first order phase transition seems to rule out the possibility of cos¬ 
mological strange quark matter. But strange quark matter, if stable, is almost unavoidable in 
dense stellar objects. Pulsars and compact x-ray sources will contain sttange stars containing 
bulk sttange quark matter ^T] ], and a number of observational signatures have been 

suggested to distinguish these from compact stars made of ordinary hadronic matter. This 
is an interesting story in its own right (see ^ for reviews), but the focus of the present 
presentation is the possibility of directly observing debris from collisions of sttange stars in 
the form of strangelets in cosmic rays. 


In the following I will summarize the main properties of strangelets, describe the propagation 
of strangelets in cosmic rays, estimate the flux to be expected in cosmic ray detectors, and discuss 
the possible detection in a couple of future experiments. Finally I will argue that strangelets may 
also be relevant for the puzzle of ultra-high energy cosmic rays. 


2. Strangelet properties 


Neglecting quark Cooper pairing (to which we return in a moment), and approximating the up 


and down quark current masses to zero, the mass of strangelets in the MIT bag model |22] depends 
on three model parameters, namely the sttange quark mass, the bag constant, B, and the strong 
fine-structure constant, o;^. In most treatments, the strong coupling and therefore is set equal to 
zero. This is clearly not correct at the relevant densities, but it has been shown [Q], that a nonzero 
value can be mimicked very well by a scaling of the bag constant, and taking = 0 simplifies 
calculations. A value of B between = 145MeV and = 165MeV allows bulk sttange quark 

matter to be absolutely stable for not-too high sttange quark mass. In fact, for < 145MeV even 
2-flavor quark matter becomes absolutely stable, having energy per baryon lower than 930MeV, 
which is ruled out by the stability of nuclei. Such a transformation would not require any weak 
conversion of quark flavor, whereas a similar decay into strangelets requires almost one-third of the 
quarks transformed simultaneously into sttange quarks, a highly suppressed process. 

Masses of strangelets can be calculated in the MIT bag model by explicitly solving the Dirac 
equation with bag model boundary conditions, filling up, down, and strange quark energy levels in 
the way which minimizes the mass for a given A. Apart from details of closed shells the outcome of 
such results can be reproduced fully in a multiple reflection expansion including volume, surface, 
and curvature terms 0, ^ 2^ . 

For a sttange quark mass of nts = ISOMeV and a bag constant B*/^ = 145MeV, the mass of a 


strangelet in the lowest energy ground state as a function of baryon number is given by |26, 271 


M(A) = 874MeVA -f 77MeVA2/^ -f 232MeVA^/^ 


(2.1) 


Thus, for these parameters bulk quark matter is bound by (930 — 874)MeV = 56MeV per 
baryon. Notice that the mass is composed of a bulk term proportional to A, a surface term pro¬ 
portional to surface area or A^/^ , and a non-negligible curvature term proportional to radius or 
A1/3_ For this specific choice of parameters, only strangelets with baryon number A > 23 are stable 
relative to nuclei. Increasing the strange quark mass and/or the bag constant moves the stability cut 
to higher A. Compared to mass formulae for nuclei, the most striking features are the presence of 
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a significant curvature energy and the lack of a significant Coulomb energy. Lack of the latter also 
means the absence of a minimum in M/A as a function of A. 

The lack of a significant Coulomb energy is due to the fortuitous cancellation of charge +2c/3 
up quarks and charge —c/3 down and strange quarks in strange quark matter with equal numbers 
of the three quark flavors. Because of the non-zero s-quark mass the cancellation is not perfect. 
Typical strangelets have slightly fewer strange quarks compared to up and down, and therefore the 
net charge is slightly positive. A typical model result (to be compared to Z Ri 0.5A for nuclei) is 


= 0.1 f 

: = 8f 


^ A 
VlSOMeV/ 
m, \ 2 


-A1/3 


( 2 . 2 ) 

(2.3) 


VlSOMeV/ 

for A -C 700 and A ^ 700 respectively (the slower growth for higher A is a consequence of charge 
screening). 

Thus a unique experimental signature of strangelets is an unusually high mass-to-charge ratio 
compared to nuclei. 

In recent years it has been shown, that quark matter at asymptotically high density has an 
interesting property called color superconductivity 3C]. Even the weakest attraction (and such 
attractions exist in QCD) leads quarks of different colors and flavors to form pairs, much like 
Cooper pairs in a superconductor, except that the binding in QCD is caused by a direct attraction 
channel rather than via indirect phonon interaction. The binding energy of a pair can be very large, 
ranging from MeV to over lOOMeV. In general these systems are called color superconductors. If 
all colors and flavors pair in an equal manner one talks about color-flavor locking. 

While the phenomenon of color-flavor locking seems generic in the infinite density limit, the 
properties of strange quark matter at densities of order or somewhat higher than nuclear matter 
density is at the focus of much current research and discussion. This is the density regime of 
relevance for strangelets, strange stars, and for quark matter cores in hybrid stars (the analogs of 
neutron stars if quark matter is metastable so that it forms above a certain density in compact star 
interiors). An additional binding energy per baryon of (10— 100)MeV is not unrealistic in these 
systems, thus significantly increasing the likelihood of absolutely stable strange quark matter and 
strangelets. 

Since Cooper pairing involves quarks with equal (but opposite) momenta, the natural ground 
state of a color-flavor locked system has equal Fermi momenta for up, down, and strange quarks. 
Equal Fermi momenta implies equal number densities, and therefore a total net quark charge of 


zero for a bulk system piQ. However, a finite size strangelet has a surface suppression of massive 
strange quarks relative to the almost massless ups and downs (massive particle wave functions are 
suppressed at a surface), so that the total charge of a color-flavor locked strangelet is positive and 
proportional to area [^] : 


Z = 0.3 


m. 


A^l\ 


(2.4) 


150MeV. 

This phenomenon persists even for very large bags, such as strange stars, so color-flavor locked 
strange stars would also have a positive quark charge [Q, 33, 34]. 
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3. Astrophysical strangelet production 


Strangelets may be produced when two strange stars in a binary system collide due to loss 
of orbital energy in the form of gravitational radiation. A strange star-black hole collision may 
also release lumps of quark matter. If strange quark matter is absolutely stable all compact stars 
are likely to be strange stars [35, 36], and therefore the galactic coalescence rate will be the one 
for double neutron star binaries recently updated in [ ^ ] based on available observations of binary 
pulsars to be 83.0^gg^j^ Myr^^ at a 95% confidence interval, thus of order one collision in our 
Galaxy every 3,000-60,000 years (but see [^8|] for a somewhat lower rate estimate). 


Each of these events involve a phase of tidal disruption of the stars as they approach each other 
before the final collision. During Ibis sfage small fractions of fhe fofal mass may be released from 
fhe binary sysfem in fhe form of sfrange quark mailer. No realislic simulation of such a collision 
involving Iwo sfrange slars has been performed lo dale. Newtonian and semirelafivislic simulations 
of fhe inspiral of sfrange slars and black holes do exisl bul fhe physics is loo differenl 

from fhe sfrange slar-slrange slar collision to be of much guidance. Simulations of binary neulron 
slar collisions, depending on orbifal and olher paramelers, lead lo fhe release of anywhere from 
10^^ — IO^^Mq, where Mq denotes fhe solar mass, corresponding lo a lolal mass release in fhe 
Galaxy of anywhere from 10^^° — 3 x per year wilh fhe collision rale above. Given fhe 

high sliffness of fhe equalion of slate for sfrange quark mailer, sfrange slar-slrange slar collisions 
should probably be expected to lie in Ihe low end of fhe mass release range, so Ihe canonical inpul 
for Ihe following calculations is a galactic production rate of 


M= lO^^Voyr^^ 


(3.1) 


All slrangelels released are assumed to have a single baryon number, A. This is clearly a huge 
oversimplification, bul Ihere is no way of calculating Ihe aclual mass speclrum to be expected. As 
demonslraled in [0 Ihe quark mailer lumps originally released by tidal forces are macroscopic in 
size; when estimated from a balance belween quark mailer surface tension and tidal forces a lypical 
fragmenl size is 

Ar^4x lo3^a2oa^o> (3-2) 

where a 2 o ~ 1 is Ihe surface tension in unils of 20 MeV/fm^ and is Ihe dislance belween Ihe slars 
in unils of 30 km. Bul subsequenl collisions will lead to fragmenlalion, and under Ihe assumption 
lhal Ihe collision energy is mainly used to compensate for Ihe exlra surface energy involved in 
making smaller slrangelels, il was shown [0 lhal a significanl fraction of Ihe mass released from 
binary sfrange slar collisions mighl ultimately be in Ihe form of slrangelels wilh A rs 10^ — 10^, 
Ihough Ihese values are slrongly parameter dependent 

The lolal flux resulls derived for cosmic ray slrangelels below are moslly such lhal values for 
some given A are valid as a lower limil for Ihe flux for a fixed lolal slrangelel mass injection if 
slrangelels are aclually dislribuled wilh baryon numbers below A. 


5 







Stable quark matter in cosmic rays ? 


Jes Madsen 


4. Strangelet flux in cosmic rays ^ 

Apart from an unusually high A/Z-ratio compared to nuclei, strangelets behave in many ways 
like ordinary cosmic ray nuclei. For example, the most likely acceleration mechanism would be 
Fermi acceleration in supernova shocks resulting in a rigidity spectrum at the source which is a 
powerlaw in rigidity as described below. Due to the high strangelet rigidity, R, at fixed velocity, 
V (R = pc jZe = Amoc^ 7 (j 8 )j 8 jZe, where p is the strangelet momentum, Amo is the strangelet rest 
mass, j3 = v/c, and 7 = (1 — strangelets are more efficiently injected into an accelerating 

shock than are nuclei with A/Z 'ta 2 (c.f. discussion of nuclei in [W]), and most strangelets passed 
by a supernova shock will take part in Fermi acceleration. 

The time scales for strangelet acceleration, energy loss, spallation and escape from the Galaxy 
are all short compared to the age of the Milky Way Galaxy. This makes it reasonable to assume 
that cosmic ray strangelets are described by a steady state distribution given as a solution to a 
propagation equation of the form 

where N{E,x,t)dE is the number density of strangelets at position x and time t with energy in the 
range [E,E + dE]. 

Given a solution for N{E ) the corresponding flux in the “average” interstellar medium with 
energies from [E,E + dE] is given by 


Be 

EE{E)dE = ^N{E)dE, 
4n 


(4.2) 


and the corresponding flux in terms of rigidity is 


ER{R)dR = ZepEE{E{R))dR 


(4.3) 


(using dE/dR = Zcj3). 

Like other charged cosmic ray particles strangelets are influenced by fhe solar wind when 
enfering fhe inner parts of fhe Solar System. The defailed inferaefions are complicafed, buf as 


demonsfrafed for nuclei in [|45[], a good fif fo fhe solar modulation of the cosmic ray spectrum can 
be given in terms of a potential model, where the charged particle climbs an electrostatic potential 
of order <I> = 500 MeV (the value changes by a factor of less than 2 during the 11 (22) year solar 
cycle). This effectively reduces the cosmic ray energy by |Z|<I> relative to the value in interstellar 
space, and at the same time the flux is reduced by the relative reduction in particle momentum 
squared, so that the modulated spectrum is 


^mod(£') = 


R{E) 




Er(E + \Z\<4>) 


(4.4) 


Solar modulation significantly suppresses the flux of charged cosmic rays af energies below a few 
GeV and effectively works like a smoofh cufoff in flux below kinetic energy of order |Z|<I>. Since 

* Section 4 closely follows Ref. |p3|, but the treatment has been slightly simplified by excluding less important terms 
in the propagation equation. 
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strangelets have a high mass-to-charge ratio they are nonrelativistic at these energies, which corre¬ 
spond to rigidities of /?sm ~ where <I> 5 oo = (500MeV). 

For cosmic rays to reach the Earth or an Earth-orbiting detector like the Alpha Magnetic Spec¬ 
trometer on the International Space Station, the rigidities have to exceed the geomagnetic cutoff 
rigidity, which is a function of detector position, and for an orbiting observatory like AMS the value 
varies from 1-15 GV as a function of time. The geomagnetic cutoff rigidity for low mass strangelets 
is comparable to or higher than the solar modulation cutoff, whereas high mass strangelets experi¬ 
ence solar modulation already at rigidities above the geomagnetic cutoff. 

Eor a non-magnetic body like the Moon, there is no corresponding cutoff, and the total flux is 
given by Emod- This makes the lunar surface an interesting laboratory for a strangelet search (c.f. 
Section 5.2; see also Ref. p 6 []). 

Given the significant uncertainty in input parameters a simple but physically transparent model 
for strangelet propagation was chosen in [|43|]. ^ is given by the following sum of a source term 
from supernova acceleration, a diffusion term, loss terms due to escape from the Galaxy, energy 
loss, and spallation. 


dt 


dN 


, + DW^N + 


dN 


5 dN 

escape T [d (A jAJ H | spallation • 


(4.5) 


The individual terms will be defined and discussed in the following. Eurther terms describing, e.g. 
decay and reacceleration of strangelets due to passage of new supernova shock waves, are discussed 


in [|43|], but will be neglected here for simplicity. 

The strangelet spectrum after acceleration in supernova shocks is assumed to be a standard 
powerlaw in rigidity with index —2.2 as derived from observations of ordinary cosmic rays. The 
minimal rigidity is assumed to be given by the speed of a typical supernova shock wave, fisN ~ 
0.005, so R mi n = y{^sN)^SNd^mQC^/Ze 5MVA/Z. The maximal rigidity from acceleration in 
supernova shocks, Rmaxj is of order 10^ GV, but the actual number is irrelevant since Rmax R mi n 
and the rigidity spectrum steeply declining. Eor a total production rate of M = lO^^°M 0 yr^' of 
baryon number A strangelets spread evenly in an effective galactic volume V, the total source term 


IS 


G{R) = 1.2 


M 


VAtuqR^ 


R 

Rmii 


- 2.2 


or in terms of energy (using dE/dR = Ze[5 and G{R)dR = G{E)dE) 

dt ^ ’ Ze/3 ■ 


(4.6) 


(4.7) 


The terms DN'^N + escape ^ where D is the diffusion coefficient, describe cosmic ray diffu¬ 
sion in real space and eventual escape from the confining magnetic field of the Galaxy. Charged 
cosmic rays are spiralling along field lines in the weak galactic magnetic field, but due to the 
very irregular structure of the field, the particles scatter on magnetic “impurities”, and the motion 
is best described in terms of diffusion. Erom studies of nuclei it is known, that cosmic rays are 
confined to move in a region significantly larger than the galactic disk, where most stars and inter¬ 
stellar matter is located. A typical value for the effective galactic volume confining cosmic rays is 
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V = 1000 kpc^. The standard leaky box approximation assumes D = 0 and ^ [escape = /"fescape) 

where Tsscape(^,Z,E) is the average escape time from an otherwise homogeneous distribution in 
the galactic volume, V. From studies of cosmic ray nuclei the escape time is known as a function 
of rigidity, R , as 


‘'escape 


8.09 X lO'^y 
nfi 



(4.8) 


where Rq = 4.7GV, 5 = 0.8 for R < Rq, and 5 = —0.6 for R > Rq. n denotes the average hydrogen 
number density per cubic centimeter {n Ri 0.5 when averaging over denser regions in the galactic 
plane and dilute regions in the magnetic halo). 

The term in the propagation equation ■^[b{E)N], describes the influence of energy loss pro¬ 
cesses. The energy loss rate b{E) = —dE/dt can be treated as a sum of ionization losses (from 
interaction with neutral hydrogen atoms and molecules). Coulomb losses (from interaction with 
ionized hydrogen), and pion production losses from inelastic collisions at high relativistic y-factor 
(threshold at 7 = 1.3). The various contributions are described in [^]. At speeds close to the speed 
of light the ionization loss is simply proportional to nZ^. 

Like nuclei strangelets have a roughly geometrical cross section proportional to for spal¬ 
lation in collisions with interstellar matter (mainly hydrogen). The corresponding spallation time 
scale is taken to be 

2 X lO^y 2/3 

nj8 


"^spallation — 


(4.9) 


At low kinetic energy the cross section can vary somewhat due to resonances etc, but such com¬ 
plications will be neglected here, since the detailed physics is unknown in the case of strangelets. 
We have also neglected the slight reduction in geometrical area of strangelets relative to nuclei due 
to their slightly larger density. The largest uncertainty in the treatment of spallation is the fact that 
strangelets (like nuclei) are not always completely destroyed in a spallation reaction. In addition to 
nucleons and nuclei smaller strangelets may result from this type of reaction, but we are ignorant 
of the physics to an extent where it is impossible to include this effective feed-down to lower A in 
a meaningful manner. Therefore, spallation is assumed to be a process destroying strangelets, i.e. 


T I spallation — • (4.10) 

ot ’t'spallation 

This leads to an overall underestimate of the strangelet flux. 

Strangelet energies are redistributed according to the propagation equation. Some leave the 
Galaxy or are destroyed by spallation. Occasionally strangelets get a new kick from a passing 
supernova shock, and in a first approximation they regain the source term relative distribution of 
rigidity. The time scale between supernova shock waves passing a given position in interstellar 
space is of order %sn ~ 10^ Y- This scale is comparable to or larger than the time scales for energy 
loss, spallation, and escape from the Galaxy, so reacceleration of cosmic ray strangelets has only a 
moderate influence on the steady state distribution. By adding energy (on average) to the particles 
it actually reduces the total flux of strangelets somewhat because higher energies make destruction 
and escape more likely (see Ref. ||^ for details). 
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Introducing the terms discussed above the steady state equation dN/dt = 0 leads to the fol¬ 
lowing differential equation for N {E) 


b{E) 


Ie 


N{E) 

t{E,A,Z) 


G{E), 


(4.11) 


where 

\/t{E,A,Z) 

1 / ^escape “hi / "^spallation “h 1 / "^1 OSS ? (4.12) 

with 1 /Tioss = —db{E)/dE. When energy loss can be neglected, b{E) ^ 0 and |tioss| ~ In this 
limit the spectrum is simply given by 


N{E)^G{E)t{E,A,Z), 


and 1/t(£',A,Z) ~ l/'tescape “h 1 /t^spallation■ 

The general solution of the propagation equation requires numerical integration, but several 
limits can be treated analytically and provide a physical understanding of the full numerical solu¬ 


tions found in [|43|]. The special cases (disregarding solar modulation and geomagnetic cutoff) can 
be divided according to the relative importance of the different time scales, Tescape, 'fspallation^ and 
I'^iossl (the energy loss time scale is negative at low and high energies, describing a net increase in 
number of particles). 

When one of the time scales is significantly smaller than the others at a given energy (rigidity), 
the corresponding process dominates the physics. The relative importance of the processes depends 
on strangelet properties A,Z, on the density of interstellar hydrogen n (though most processes have 
the same n-dependence), and of course on the strangelet energy, E (or rigidity, R). Energy loss 
dominates at low energy, spallation at intermediate E, and escape from the Galaxy at the highest 
energies. 

The energy loss domination is important only at energies below the solar modulation cutoff. 


so I refer the reader to |43] for details. 


At intermediate energies the spectrum is determined by the strangelet spallation time, T 
Tspallation^ and the energy distribution is approximately given by 


N{E) ^ G(E)Tspallation(«,j8,A) 

with Tspaiiation = 2 X lO’yn^' 

This gives the approximate result 


(4.13) 


T)?(/?) = 2.34 X lO^m ^yr ^sterad" 




(4.14) 


or for the total flux above rigidity R 

E{> R) = 1.95 X 10^m^V^^sterad^'A^°'^'^’^Z^'-2/?Q^-^A. 
The results scale in proportion to 

jSsN \ / 0.5cm~^ \ I M 

0.005 ) \ n J \ lO-iOMeyr-i 



/lOOOkpc^X /930MeV\ 
V V y V moc^ J 


(4.15) 


(4.16) 
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Notice that the differential strangelet spectrum keeps the source term slope, G{R) oc 

The intermediate energy domain is replaced by the high energy domain when Tescape < '^spallation- 
Except for very low A this happens when R > l.OGVA' or E > l.OGeVZA' At high en¬ 

ergies the spectrum is determined by the confinement time of strangelets in the Galaxy, T 'Tescape. 
and the propagation equation leads to 

N{E) G(E)Tescape(«,i8,7?). (4.17) 

For semirelativistic or relativistic strangelets with j8 Ri 1, Tescape so the spectrum is steep¬ 

ened from the source term R ^ ^ to 
This gives the approximate result 

Fr{R) = 2.40 X lO^m^V^'sterad^' (4.18) 

and for the total flux above rigidity R 

F{> R) = 1.33 X 10^m^2yj.-istgj.ad-i^o.22-i.2^^L8^ ^4 ^9^ 

The astrophysical input parameters in the present calculations are uncertain at the order of 
magnitude level, so approximate relations for the total flux of strangelets hitting the Earth or Moon 
accurate within a factor of 2 (for fixed input parameters) are useful. As indicated above solar 
modulation effectively cuts off the strangelet flux at rigidities of order /?sm ~ 
which is in the part of the spectrum where the strangelet flux is governed by spallation. The total 
flux hitting the Moon or Earth is therefore roughly given by 

Etotai ~ 2 X 10^m^^yr^'sterad^'A^°"*^^Z^''^max[7?sM,^Gc] (4.20) 


depending on whether solar modulation or geomagnetic cutoff dominates. In the case of solar 
modulation domination (always relevant for the Moon, and relevant for AMS as long as /?sm > 
7?gc) one obtains 

Etotai ~ 2 X 10^m^V^‘sterad^U^' °^’Z^°-®<I>50°j*^A. (4.21) 

For strangelets obeying the CEE mass-charge relation Z = 0.3A^/^ this becomes 


Etotai ~ 4 X 10^m^2yj.-i sterad^U^*-4®'^<I>5o°o^A 

Ri 2.8 X 10^m^V^^sterad^^Z^2.2<j,-a6^^ 


(4.22) 

(4.23) 


which reproduces the numerical results of Ref. [ ]43| ] to within 20% for Z > 10 and to within a 
factor of a few even for small Z, where the assumptions of nonrelativistic strangelets and spallation 
domination both are at the limit of being valid. 

The strangelet flux for A of order unity is high enough to be of interest for various upcoming 
experimental searches, and at the same time small enough to agree with previous searches which 
have given upper limits or shown marginal evidence for signatures consistent with strangelets (see 
[ jn] ] for an overview). As stressed several times above, many parameters are uncertain at the order 
of magnitude level. The scaling with these parameters is indicated where relevant. In particular 
this is true for the overall normalization of the strangelet flux as expressed via the parameter A (Eq. 
(|4.16)), whereas the relative behavior of the differential flux is less uncertain. 
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Apart from the uncertainty in parameters within the picture discussed here, one cannot rule 
out the possibility that some of the basic assumptions need to be modified. In addition to strangelet 
production in strange star collisions, it has been suggested that a (possibly small) flux of strangelets 
may be a direct outcome of the Type II supernova explosions [p7|], where strange stars form. The 
treatment above does not include such additional strangelet production mechanisms, and therefore 
the flux predictions are conservative. Another assumption that leads to a conservative lower limit 
on the flux is that spallation is assumed to destroy strangelets completely. At least for low-energy 
collisions one would expect that fragments of strangelets would survive, but the input physics 
for performing a realistic strangelet spallation study is not sufficiently well known, and there¬ 
fore the conservative assumption of complete destruction was made. A numerical simulation of 
strangelet propagation in Ref. p8[ ] assumed stripping of nucleons rather than complete destruction 
of strangelets in interstellar collisions and studied two specific sefs of values for fhe mass and en¬ 
ergy of sfrangelefs injected info fhe interstellar medium. Several ofher assumptions made in fhaf 
numerical sfudy differ from fhose of fhe currenf invesfigafion, so a direcf comparison of fhe re- 
sulfs is nol possible, excepf fhaf bofh sfudies are consisfenf wifh fhe possibilify of a significanf, 
measurable sfrangelef flux in our pari of fhe Galaxy. 

Buf ulfimafely fhe question of whefher sfrangelefs exisl in cosmic rays is an experimenlal issue. 

5. Experiments underway 

Several experimenfs have searched for sfrangelefs in cosmic rays and/or have had fheir dala 
reanalyzed fo provide limifs on fhe sfrangelef abundance. While some interesting evenfs have been 
found fhaf are consisfenf wifh fhe predicfions for sfrangelefs, none of fhese have been claimed as 
real discoveries. Interpreted as flux limifs ralher lhan detections fhese resulls are consisfenf wifh 
fhe flux eslimales given above. For discussions see iHl- 

If fhe interesting evenfs were aclual measuremenls, Iwo experimenfs fhaf are currenlly under¬ 
way (wifh fhe aufhor as a humble fheorefician participating in collaborafions wifh very capable 
experimenfalisfs!) will reach sensitivities, fhaf would provide real sfafisfics. These experimenfs are 
fhe Alpha Magnetic Specfromefer (AMS-02) on fhe Infernafional Space Sfafion, and fhe Lunar Soil 
Sfrangelef Search (LSSS) al fhe Wrighl Nuclear Sfrucfure Laboratory af Yale. 

5.1 AMS-02 

The Alpha Magnefic Specfromefer (AMS) is a high profile space-based parficle physics exper- 
imenf involving approximately 500 physicisfs from more lhan 50 insfilufions in 16 counfries, led 
by Nobel laureate Samuel Ting of fhe Massachuseffs Insfifufe of Technology (MIT). The AMS flew 
in space in June of 1998 aboard fhe Space Shuffle Discovery and if is currenfly scheduled fo 
fly to fhe International Space Sfafion (ISS) on Utilization Flighf #4.1 (UF4.1) wifh launch in 2008. 
Once on fhe ISS, fhe AMS-02 will remain active for al leasl Ihree years before if is relumed fo 
Earlh. AMS-02 will provide dala wifh unprecedenled accuracy on cosmic ray eleclrons, positions, 
protons, nuclei, anti-nuclei and gammas in fhe GV-TV range in order fo probe issues such as an- 
limaller, dark mailer, cosmic ray formation and propagafion. And in addition if will be uniquely 
suifed fo discover sfrangelefs characterized by extieme rigidilies for a given velocily compared fo 
nuclei O- AMS-02 will have excellenf charge resolution up fo Z Ri 30, and should fherefore be 
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able to probe a large mass range for strangelets. A reanalysis of data from the 1998 AMS-01 mis¬ 


sion has given hints of some interesting events, such as one with Z = 2,A = 16 []50|], but with the 
larger AMS-02 detector running for 3 years or more, real statistics is achievable. 


5.2 LSSS 

The Lunar Soil Strangelet Search (LSSS) is a search for Z = 8 strangelets using the tandem 
accelerator at the Wright Nuclear Structure Laboratory at Yale. The experiment involves a dozen 
people from Yale, MIT, and Arhus, led by Jack Sandweiss of Yale. The experiment which is 
currently in its calibration phase, studies a sample of 15 grams of lunar soil from Apollo 11 (sample 
number AH10084). The expectation is to reach a sensitivity of 10^^^ over a wide mass range. 


sufficient to provide detection according to Equation (4.21) if strangelets have been trapped in the 
lunar surface layer, which has an effective cosmic ray exposure time of around 500 million years. 


6. Moving the GZK-cutoff-Strangelets as ultra-high energy cosmic rays 


One of the most interesting puzzles in cosmic ray physics is the apparent existence of cosmic 
rays with energies well beyond lO^^eV, with measured energies as high as 3 x lO^^eV. Briefly 
speaking there are two puzzles |M, ^]: 


1. It is almost impossible to find an astrophysical mechanism capable of accelerating cosmic 
rays to these energies. 

2. Even if acceleration happens, ultra-high energy cosmic rays have a relatively short mean-free 
path for interactions with the cosmological microwave background photons, and only cosmic 
rays from fairly nearby (unknown) sources would be able to reach us with the high energies 
measured. 


Interestingly, strangelets circumvent both of these problems, and therefore provide a possible 
mechanism for cosmic rays beyond the socalled Greisen-Zatsepin-Kuzmin (GZK) cutoff. The 
benefits of strangelets with respect to the two problems are []5^ : 


1. All astrophysical cosmic ray “accelerators” involve electromagnetic fields, and the maximal 
energy that can be transmitted to a charged particle is proportional to the charge, Z. Eor 
instance some mechanisms can provide a maximal rigidity, /?max> proportional to the product 
of the magnetic field and size of the region where acceleration takes place, and for relativistic 
particles this corresponds to a maximal energy Emax = -ZEmax- The charge of massive 
strangelets has no upper bound in contrast to nuclei, so highly charged strangelets are capable 
of reaching energies much higher than those of cosmic ray protons or nuclei using the same 
“accelerator” |53]. 


2. The GZK-cutoff is a consequence of ultrarelativistic cosmic ray projectiles hitting a 2.7K 
background photon with a Eorentz-factor y large enough to boost the 7 x lO^^eV photon to 
energies beyond the threshold of a process leading to significant projectile energy loss, such 
as photo-pion production, photo-pair production, or photo-disintegration. The threshold for 
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such a process has a fixed energy, f’xhreshoid. in the frame of the cosmic ray (e.g., f’xhreshoid ~ 
lOMeV for photo-disintegration of a nucleus or a strangelet), corresponding to a Lorentz- 
factor Trhreshoid = £^xhreshoid/£^ 2 . 7 K (~ 10^° for photo-disintegration). But this corresponds to 
a cosmic ray total energy 

^Xotal T'rhreshold^nJoC . (6-1) 


Thus, f'xotai (~ lO^^A eV at the photo-disintegration threshold) is proportional to the rest 
mass of the cosmic ray, and since strangelets can have much higher A-values than nuclei, 
this pushes the GZK-cutoff energy to values well beyond the current observational limits for 
ultra-high energy cosmic rays [^, 541. 


A testable prediction of the strangelet scenario for ultra-high energy cosmic rays is that strangelets 
at a given energy will be more isotropically distributed than protons and nuclei because the gyro- 
radius in the galactic magnetic field is proportional to so the arrival direction for high-Z 
strangelets points back to the source to a lesser extent than for low-Z candidates which have gyro- 
radii comparable to the size of the Galaxy. 


7. Conclusion 

The total strangelet flux reaching the Moon or a detector in Earth orbit is in a regime that could 
be within experimental reach, and therefore provide a crucial test of the hypothesis of absolutely 
stable strange quark matter. Experiments are underway which are sensitive to the high mass-to- 
charge signature expected for such events within the tlux-range predicted. Strangelets may also 
provide an interesting explanation of ultra-high energy cosmic ray events. 
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